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ABSTRACT: 



The report presents a computer program for prediction of performance 
of single-stage axial tixrbines of given geometry. The three-dimensional 
method developed by Vavra is applied, taJcing account of streamline 
curvatures and slopes, as well as enthalpy and entropy gradients in 
the solutions of the equation of motion, and of boundary layer thick- 
nesses in the continuity equation. 

A choice among five different loss correlation methods and two 
flow angles correlations is offered. Loss coefficients and flow 
angles are automatically calculated from blading geometry and actual 
flow conditions for every streamline, according to the selected 
correlation method. 

A fair agreement of predictions with several actual turbines 
experimental results was found in ref, [4], where also the applicability 
of different available correlations is discussed. 
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SUBSCRIPTS 



A Axial 

E Equivalent 

is Isentropic expansion from total inlet conditions 

m Mean streamline 

m Meridional 

R Rotor 

r Radial direction, cylindrical coordinates 

S Stator 

z Axial direction, cylindrical coordinates 

9 Peripheral direction, cylindrical coordinates 

0 Station ahead the stator 

1 Station between stator and rotor 

2 Station after the rotor 
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INTRODUCTION 



An accurate knowledge of performance achievable with turbines 
of fixed geometry, both at design and off-design conditions, is 
important in advanced propvilsion and power systems. The more accurate 
the prevision is, the more meaningful the optimization of the txarbine 
itself and of the overall system will be . 

A first way to solve this problem consists in the testing of a 
prototype, but it presents several disadvantages, namely it is 
es^ensive and time consuming. Furthemore, if effects of turbine 
geometry changes have to be investigated, different prototj^es 
will be required. 

An alternative solution to the problem consists in theoretical 
methods of prediction. The present availability of high speed 
computers allows one to solve systems of differential equations and 
to perform a large number of iterations in a short time. Moreover 
computer programs give the possibility of changing in^iut data, so that 
a large variety of different turbine configurations and of working 
conditions can be investigated. 

This report presents a computer program, derived by the 
three-dimensional calculation method given first by Vavra [l] and 
then developed by Eckert C 2 ] and Harrison This method takes 

account of enthalpy and entropy gradients and of streamline curvatures 
and slopes in equations of motion, and of boundary layer thicknesses 
in continuity equation. It is applicable to single-stage axial 
turbines with unchocked flow. 

A critical point in all theoretical performance calculations 
is the prediction of the losses that occur in the blade rows. In 
the present program, one can choose among five different loss correla- 
tion methods. For each method pertinent subroutines compute loss 
coefficients from blading geometry and flow characteristics. A 
comparative study of these methods can be found in [4], where they 
were applied to actual turbines of different design practices, in a 
wide range of all significant flow and geometrical properties. 
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The comparisons presented in [ 4 ] with experimental performance 

showed a fair agreement of predictions with all turbine stages 
that were examined. 

The present report is divided into four sections: in the first 
the used equations are set up and necessary assumptions are discussed; 
in the second the method of solution is described; in the third 
different losses and discharge flow angles available prediction methods 
are described; and in the fourth the instructions for the use of 
conputer program are given, with an example of imput and output data. 
Eventually in the appendix the FORTRAN list of the MAIN routine and 
of the 28 subroutines is given. 

METHOD OF ANALYSIS 

The flow in tirrbomachines is ruled by the three fundamental 
laws of fluid dynamics, namely: 

Equation of Motion 

Energy Equation 

Equation of Continuity. 

All proper three-dimensional calculation methods proposed 
in the technical literature solve the differential equations derived 
by these laws. Differences among these methods are due to the three 
following reasons: 

different stations are chosen where equations are solved; 
different hypotheses are formulated to describe the real 
processes occurring in the machine; 

different numerical methods are applied in the solutions of 
the equations. 

As far as the first problem is concerned in the present method, 
the equation of motion is applied at three stations, namely ahead of 
the stator, after the stator, and after the rotor, indicated respectively 
with subscripts 0, 1, 2 in Fig. 1. On the contrary, since blade geometry 
is given, the continuity equation is applied at blade throat sections 
1* and 2* in Fig. 1, as suggested in [ 5 ] . 
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those occurring in the boundary layers along the blade surfaces, and 
flow frictional losses between blades can therefore be neglected. 

Assumption k) is reasonable for uacooled blades, where energy 
changes due to heat transfer are negligible with respect to those 
related to velocity changes. If blades are cooled, the additional 
term q in the energy equation must be specified. Assumption 5) is 
valid in most cases. If the thermodynamic processes occur in ranges 
where real gas effects are not too strong, average values of gas 
properties give good approximations. 

Some additional assumptions, regarding streamlines slope and 
curvature, are discussed in the section pertaining to the equation 
of motion. 

As far as the third problem is concerned, the numerical 
small difference method was found to give a good convergence, and 
was therefore applied for the solution of the differential equations. 

The following sections give the derivations of the equations 
with the above-listed assumptions. 

Equation of motion 

The general equation of motion for relative flow in vectorial form, 
is: 



St? . 

•^+VHr = Wx(VxW+2o)) + TVS + 



( 1 ) 





j 0) = 0 


(2) 


For absolute flow, where ' 


( V = W 


(3) , equation (l) 




II 


(4) 



becomes: 



|X+VH=Vx (VxV) + TVS+ff 
With hypothesis 
and hypothesis 






2 ) = 0 



(la) 
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This allows one to more accurately evaluate the pressure drops 
occurring in blade channels, since in this way it is possible to 
account for boundary layer thicknesses. This problem is particularly 
important at high Mach numbers where a small change in flow area causes 
large differences in pressure ratios. 

As far as the second problem is concerned, the following 
assumptions are made: 

1) Steady flow 

2) Frictionless flow at the stations where the equations of 
motion are solved. 

3) Axisymmetric flow at the stations where the equations of 
motion are solved. 

4) Adiabatic flow. 

5) The working fluid is a perfect gas, namely a gas having 
the equation of state: 

p = P % T 

and specific heat at constant pressure is independent of temperature 
and pressure. 

Hypotheses l) and 3) are interrelated since the flow can be steady 
in both rotor and stator only if it is axisymmetric. Although they 
are correct only for cascades with an infinite number of blades, 
axisymmetry and steady conditions are necessary assimaptions since any 
other hypothesis would require unknown quantitative information on 
flow inside the blade rows. One problem connected with these ass;xmptions 
is that downstream conditions cannot affect the flow upstream. There- 
fore pressure distributions at rotor exit can deviate from those imposed 
by the discharge conditions. 

With assumption 2), all entropy changes are supposed to occur in 
the blade rows ahead of the considered station. This assumption can 
be justified by the fact that velocity gradients in the regions, where 
the equation of motion is applied, are in all probability smaller than 
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H I 



equations (l) and (la) become respectively: 

V Hp = W X (V X W + 2u)) + T V S (5) 

VH = Vx(VxV) + TVS (5a) 



The following derivations for equation (5) are valid also for 
equation (5a), if substitutions (2), (3), (4) are performed. 

It is convenient to express equation (5) in cylindrical coordinates 
(see Fig. 1 for symbols). Equating the three components of equation (5), 
the following equations are obtained: 



a. tangential component: 

dHR _ ^ _Wr I d(rWu) 

S9~r‘"d9~ cz"^~r“' r "59-^ 



T dS 

2^ r 



( 6 ) 



b . meridional component : 



dHp “ dWr SWa 1 Wu r SWa 

dz ^ dz " dr " r 59 



dz 



(7) 



c. radial component: 

Wu I d(rWu) 'SWp”! 

dr ~ r ■ dr d9 

with h3q)othesis 3) ( 

and hypothesis 4) w d t • ^ Hr = 0 

the system of equations (6), (7), (8) reduces to equation (9): 



dW dWA n . 

- "a L -Si - -57 J + 2“ »U + (8) 

) = 0 



dHR d(rW^) 

dr ~ r §7^ 



Wa 



dW^ 

dz 



dW; 



+ Wfl 



or 



+ 2uj Wu + T 



or 



(9) 
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Using substitution: 



( 10 ) 



He = Hr + ^ 



with hypotheses 5) H = Cp T and 4) Hg = const (see Fig. 2) 

TT 2 

Ht? Wp 

"2 = ^ - -iz 



( 11 ) 



and introducing the non-dimensional quantities: 



Wfl 



Y = 



W 



Am 



(12) 



X = 

^m 



(13) 



S* = |- 



(14) 



Equation (9) becomes: 

1 dY^ ^ 2 _ . 2rjn dW^ p 

COS p ( - -r; ^ “ 

y2 ox Oz 



dS*N 



dp 



o — ^ ) + 2 tan P 

2 cos2x ' P dx 



2 . 2„ 4 sin p cos p ^ Uj^ U cos p 2 cos^ P 



+ — sin P + 

X 



^Am Y 



Wa^ Y^ 



WAm^ ^ 



CX 



r2 Hp cos2p 

L-Vr- 

V ^ 



- Sin 



in^ J 



dS* 

dx 



= 0 



( 15 ) 



In this equationj terms containing cos A and ^ — depend respectively 

o z 

on streamlines slope and curvature. Vavra [iJ showed that the most 
reasonable assumptions, when only three points of a streamline are known, 
are the following (see Fig. 1 for symbols): 



Wa dz - j2 



(16) 



p 

cos X = 



l 2 t (§5)2 



l4 



(17) 



In equation (l6) the positive sign holds for station (2) and the negative 
one for station (l), for 6r positive as shown in Fig. 1. Suggested values 
for k are between 4 and 6. Finally with equations (l6) and (17) into (15) 
the following equation is obtained: 



ARs2 



(in y2) r l2 + (|«) 

= - cos^p L 2 k r„ 

dx ^ - m £2 

d6 2 2 ^ sin3 cosP 

- 2 tan 6 - sin p - 

■ dx X 



dS* 

dx 



Wm y 



2 U cos^p 






2 cos^p dHg 



-2 Hg cos^ „ 

L - — ■ > ” - P 

^Am ^ 



dS-^ 

dx 



(18) 



For the station after the stator row, the following equation can be 
derived with (2), (3) and (4) in equation (I8): 



d(ln Y^) _ 
dx 



= - COS 



O' 



(.2 K 






o / 



(AR/2)^ 



dS^ 

dx 



-2 tan 0^2^ 



dcy 

dx 



2 . 2 

— sin 

X 



di + 



2 cos‘ 



O' 



y 2 V„ 2 
Am 



M j 
dx " ^ 



2 H cos2 a 
y 2 V^2 



. 2 

sin 




dS* 

dx 



(18a) 



Equation (18a) could also have been derived directly from equation (5a). 
In equations (I8) and (l8a) the last term is depending on flow 
irreversibilities due to frictional losses. The specific entropy S* will 
be calculated, after Cl], by: 






S* = In 



2 Kg cos"^p 



^ Y 






2 Hg cos p (1-;) 



(19) 
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and by: 



2 2 



Y V 



Am 



S* = In 



2 Hq cos2p 

2 Ho cos2p(l-C) ^ 



+ S, 



(19a) 



where Q represents the kinetic energy loss coefficient which is 
established by the chosen loss correlation method. 

Energy equation 

The energy equation in a relative flow is given by: 



^dt'V (Hr) = l?dt''7(q) 



a 



(■^) dt 



( 20 ) 



and in an absolute flow, with (2), (3) and (4) in (20), : 
7dt-v (Hr) = ’^dfv(q) . ^ (.L) dt 



(20a) 



where q is the heat added to a unit mass particle from external 
sources. With assumptions l) d/^( ) = 0 and 4) q = 0, equations 
(20) and (20a) become: 

V?dt-V (Hr) = 0 
^dr-V (H) = 0 

Equations (2l) and (21a) indicate that for steady, adiabatic flow 
the relative total enthalpy is constant along any streamline. This 
does not mean that V (Hr) = 0, since relative total enthalpy can assimie 
different values for different streamlines. 

Equation of continuity 

Pressure ratios across a stage necessary to discharge a given 
flow rate are function of inlet Mach number, blade openings and the 
losses occurring in the blades. For an isentropic flow the mass flow 
rate passing thru a section dAj^ perpendicular to the flow is given by: 



( 21 ) 

(21a) 



l6 



(22) 



d 



m 



IS 



^is ^is 



where P^g and are respectively the density and the velocity 
corresponding to an isentropic process. It can be shown that 
equation (22) can be expressed for a perfect gas (see assumption 5) 
by: 

d-J^iis = $ dAn (23) 

v/R(j Tip 



where Pp, Tp are, respectively, total pressure and total temperatixre 
(constant for isoentropic and adiabatic flows) and f is a function of 
pressure ratio and gas isentropic exponent, namely; 



$ = 





2/y p Y+11 

i"p 



(24) 



With the symbols in Fig. 3, assuming the pressure to be constant along 
the throat "a" equation ( 23 ) becomes: 

d liiis = ^ ^ ^ d r ( 25 ) 

/ Rq Tp 



Integrating equation ( 25 ) along the radius and for z blades the 
flow rate passing thru the stage is: 



“is 



Rp 

" J 
% 



% Tp 



# a cos A. dr 



( 26 ) 



In actuality, the flow rate will be less than %g, because of the losses 
occurring in the blade channels. 

A so-called restriction factor; 

p dm 
d mis 

is now introduced. 
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In [ 5 ] the following relationship between § and the kinetic 
energy loss coefficient is derived: 






1 + 

“ K 



(28) 



where H and K are respectively the boundary layer wake form factor 
and the energy factor. For boundary layer profiles following the 
power law, the ratio -jr can be expressed in power series: 



H 

K 



1 - (l-Xe) ^ LXe^"Vl+(2i-l)m 



(29) 



(1-Xe) Zi |_Xe^'Vl+(2i-l)m J- (l-Xe) Lxe^'Vl+(2i+l)m 



where : 

Y-1/ Y 

Xe = 1 - (P/Py) ( 30 ) 

and m is the exponent of the power law for the boundary layer profile 
(for incompressible flow, m depends on Reynolds number and varies between 
1/7 and 1/10). Experience has shown that the loss coefficient C* has 
to be taken less than the value C^OT that corresponds to the overall 
energy losses that are given by the correlation methods. In the present 
program one of the three following different assumptions for C* can be 
chosen: l) C* = CrpQrp ; 2) C* = -ICtot ’ 3) C* = CpROFILE • 

While assumption l) is certainly conservative, it is difficult to state 
which of assumptions 2) or 3) is the best one. Good agreement between 
experimental values and predicted flow rates for all considered turbines 
at low tip clearance values was found in [ 4 ^, where assumption 2) was 
adopted. Assumption 3) was however foiuid in closer agreement with 
test data for large values of rotor tip clearances. 
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Tip leakage flow rate will be calculated by using the relation: 



m, 



CL 






V R Trr 



cos Acl 



(31) 



where: = 2tt c^ 



(32) 



Since the flow is not perpendicular to the leakage area, and since the 
restriction factor C,pjp does not account for boundary layer blockage 
occurring on shroud, the leakage flow rate must be divided by a 
coefficient (larger than l)^ assumed equal to 2 in the present 
program. 

With equations (29) and (33), the overall flow rate passing thru 
a blade row is given by: 



“tot = “ + “cl = ^ 



Rt 

J 



R 






/ 



• rp 2. 

$ a cos X dr + — ( § 

A, 



(33) 






CL 



/ 



R Tr, 



$ cosX)^PP A^P 



Equation (33) can be expressed in the following dimensionless form: 



A = z Si 



^ 2n » ^T c^ / T 

(?. $ a cos A ) dx + — ^>riiTP = ( ? cos A) 

^ Rg T ^CL Si za /R Tj 



(34) 



TIP 



where a = — 
a 
m 



(35) 



METHOD OF SOLUTION 



Since no direct solution of governing equations is possible, an 
iterative method is followed. At first, solutions are found without 
considering the influence of streamline curvature and slope. The 
derived radial shifts are then used for the second and final cycles of 
calculation. In Fig. 4 the scheme followed in the computer program is 
indicated. The names in brackets correspond to the names of the 
subroutines which perform the operations described next to the brackets. 
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The flow ahead of the stator is supposed to he uniform, so that the radii of 
streamlines dividing the flow rate in equal paxts can he easily computed. 

For stator and rotor rows the following functions are supposed to he 
known: 



STATOR 

(r, flow condition) 

^1 ®1 

=s “ ~s flow condition) 



ROTOR 

^2 ~ ^2 condition) 

ag = (r) 

=R - flow condition) 



All derivatives used in the following analysis are computed hy 
finite difference methods. The computer program calculates flow 
conditions along 5 streamlines. However the same scheme of calcu- 
lation can he applied to any number of streamlines. By assuming the 
inlet Mach number Mq, the flow rate passing thru section 0 can easily he 
determined. Radial position of streamlines at station 1, Rj_, and the axial 
velocity at mean streamline? are also assumed a priori, and will 

he checked later. 

For any given values of AR, and 6r, equation (l8a) can he 

solved hy iteration. In fact it can he reduced to: 



d(Xn Y^) 
dx 



I (x) 



^x = 1 

Equation (36) integrated, with boundary condition |y = ]_ 

gives: ^ 

2 r 

In Y = I dx 

tJ 

1 

and 

pX 

Y . 'Jl 



(36) 

(37) 

(38) 



(39) 
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Solutions are found by ti'iai and error, assuming approximate values 
of Y(x) and verifying equation (39)* Once equation (39) has been 
solved, a1 1 flow conditions at station (l) are known and the overall 
continuity equation (34) can be solved. If the calculated flow rate 
is different from the assumed one, will be changed and the 

iteration shown in Fig. 4 is performed. When overall continuity 
equation is satisfied, the assumed radial position of streamlines 
is checked, and necessary iterations are performed. If streamlines 
continuity is met also, the flow at the rotor inlet is known completely. 
A similar iterative method is then followed for station 2 after the 
rotor as indicated in Fig. 4. If the flow after the rotor is known 
also, the radial shift of the streamlines can be determined, and the 
whole calculating procedure is repeated by taking accoimt of the 
streamlines curvatures and the slopes in solving the equation of 
motion. 

Eventually the resiilting average pressure ratio across the turbine 
is compared with the desired one , and if the difference is larger than 
a smallprefixed error, the complete calculating scheme is repeated 
again with a new value of inlet Mach number. After the specified 
pressure ratio has been obtained, all output information is calculated 
and printed. 

FLOW ANGLES AND LOSS COEFFICIENTS CORRELATIONS 

Evidently, even the most accurate method of performance prediction 
will give results in disagreement with actuality, if the flow angles 
and/or the losses in the blade rows are not predicted with the necessary 
accuracy. 

Several correlations, either empirical or theoretical, have been 
proposed in the technical literatiire . A comparative study of five 

correlations methods available in the present program can be found in 
Reference C 4] . The purpose of the present chapter is not to discuss 
the applicability of these different correlations, but only to indicate 
the problems which arose, in applying the different correlation methods. 
In Fig. 5 the nomenclat\ire and sign convection of the velocity triangles 
is indicated. Unless otherwise stated, the variables are positive if 
directed as in Fig. 5- In Fig. 6 the blading geometry characteristics 
used in the different correlation methods are indicated. The names 



21 



given in Fig. 6 are the FORTRAN names used in the program for stator 
quantities at the mean radius. 

The program is organized in such a way that other than considered 
flow angles and loss coefficients correlations can he introduced by- 
adding a few cards and pertinent subroutines. Since the solution method 
is iterative, added subroutines can use blading and flow characteristics 
that are not used in the present program. 

a) Flow angles correlations 

Two methods, one given by Ainley and Mathieson [6J and the other 
given by Traupel [lO] are available in subroutines. In both these 
methods, exit flow angles ( ^2 )are considered to be independent 

on inlet angles ( Oq, . While in Traupel's method the angles are 
a function of blading geometry only, in A&M method also the exit Mach 
number is supposed to affect the angles. Although both these methods 
are given for a "reference diameter approach", they will be applied 
in three-dimensional calculations with the assumption that at any 
radius the angle is related only to the corresponding characteristics. 
Since during the calculations both Mach numbers and streamlines 
radii will be changed, flow angles corrections are automatically 
performed during the solution of the equation of motion. 

b) Loss coefficients correlations 

The assumption is made that at any particular radius the loss 
coefficient is a function of the characteristics at that radius only. 

It is realized that this assumption is not realistic, since secondary 
and leakage losses are certainly not equally distributed along the 
radius. On the other hand, secondary and tip leakage flows are not 
understood well enough at the present time to formulate any other 
reasonable assumption. 

Five loss correlation methods given in Ref. [6-lOj are available 
in subroutines. They calculate loss coefficients from given blade 
geometry and flow angles and velocities. 

The first considered method is the Ainley and Mathieson L6] 
method. Recently an improvement to the original method was given 
by Dunham and Came [7]. Both methods are available in subroutines. 

These methods express losses by total pressure loss coefficient Y, 
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defined by: 



P 



P 



Y = 



To 



t 



( 40 ) 



p. 



t 



- P 



and assume Y to be independent on pressure ratio Pijq/p . The 



relationship between the kinetic energy loss coefficient C and Y is 
given by: 



The assumption that Y is constant implies that C decreases with increasing 
Mach number. Since this behavior is not experimentally proved, two 
alternatives to equation (4l) are available in the program. In both 
these alternatives ^ is supposed to be constant with pressure ratio, 
and assumes respectively the value corresponding to very low Mach 
niimbers (eq. 42) or to an exit Mach number of 0.8 (eq. 43). 



Y-1 




(41) 



Y-1 





(42) 



Y-1 



^ = ^ m =.8 S 

(l/S*) Y -1 




(43) 



where p* is the pressure ratio corresponding to an isentropic Mach 
number of 0.8, given by: 



The same alternatives are given also for the simple method of 
Lenherr and Carter [9J* 

For the method of Balje ' L Oil , the boundary layer wake form 
factor H must be introduced. 

With the same nomenclature and assumptions as in eq. ( 29 ), H can 
be expressed by: 



i-1 



H = 



1 - (l-X,) 



1+m (2i-l) 



(1-Xj Ei( 



X 



i-1 



l+m(2i) 



+ (1-X^) Ei f ^e 

1 ^l+m(2i-l) 



( 45 ) 



In this program tiie assumption m = 0.15 is made. If another exponent 
has to be introduced, card #17 must be changed in the MAIN routine. 

In Traupel's method, the loss coefficients depend on the 
coefficient: 

Xp = Xp (Re, relative surface roughness, turbulence factor). (46) 

In the present program, the standard value of X =1 is assumed. If 
the considered tixrbine has peculiar values of surface roughness or 
Reynolds numbers, it will be necessary to change card #15 in the MAIN 
routine . 

For conformity with the other methods, tip clearance losses are 
expressed by the kinetic energy loss coefficient Ccl? instead of by 
the overall efficiency decrease AT|qj^ as in [ lO] . The relationship 
between ATI^l Ccl is given by: 



^CL ~ (l“^r) ^ 



LI - (1 - 



^is 2 ] 



( 47 ) 



where Ahj^g is the isentropic enthalpy drop through the tixrbine and 
is the rotor kinetic energy loss coefficient for zero tip clearance. 
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CCMPUTER PROGRAM 



The present computer program is composed of a MAIN routine 
and 28 subroutines, eight of which are concerned with flow calcula- 
tions, two with flow angle correlations, 15 with blade loss 
correlations, and three with numerical interpolations. A brief 
description of the purpose of each routine is given below. 

a) MAIN routine 

1. The primary function of the MAIN routine is to control 
the overall logic of the computer program which is indicated in 
Fig. 4. Furthemore the imput data are read, numerical constants 
are introduced, interpolations concerning stall and blade angles and 
flow areas are performed and the flow angles are computed. Eventually, 
also the output data are printed. 

b) Flow calculation subroutines 

2. In subroutine CHAW the flow rate passing through the section 
0 is calculated from the conditions ahead of the stator. 

3 . In subroutine STAT;^R the equation of motion (l8a) is solved, 
and the flow conditions after the stator are calculated. If necessary, 
new values of the flow angles are also calculated. 

4. In subroutine Rj^T/iRl the absolute flow properties calculated 
in STAT)6 r are converted into relative properties at rotor inlet. 

5 . In subroutine r/)T^R 2 the equation of motion (I 8 ) is solved, 
and the flow conditions after the stator are calculated. If necessary, 
also new values of flow angles are calculated. 

6. In subroutine FL/)¥R the flow rate passing through the blade 
sections, either at station 1 or 2, is con^juted by equation (34), and 
compared with the value calculated in CHAN. If the difference is larger 
than a small prefixed error, a new value of the axial velocity at the 
mean stremaline is computed. The fractional values of the flow rate 
passing between the hub and each streamline are calculated also. 

7 . In subroutine SLIWE the fractional values of the flow rate 
computed in FL^^WR are compared with prefixed ones. If the differences 
are larger than a small specified error, new values of the streamline 
radii, either at station 1 or 2, are computed. 
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8. In subroutine ALj6si the stator loss coefficients are 
calculated from the chosen loss correlation, as a function of 
blading geometry and flow characteristics. Coefficients C* of (28) 
are computed for the stator blades. 

9 . In subroutine AL^S2 the rotor loss coefficients are calcu- 
lated from the chosen loss correlation as function of blading 
geometry and flow characteristics. Coefficients C* of equation (28) 
are computed for the rotor blades. 

c) Flow Angles Subroutines 

10. In subroutine AITGAIN the flow angles are confuted from 
blading geometry and exit Mach numbers with the method of Ainley 
and Mathieson [6]. 

11. In subroutine ANGTRA the flow angles are computed from 
blading geometry with Traupel's method ClO]. 

d) Loss Correlation Subroutines 

The six subroutines given below are concerned with the methods C6J and 

C7J . 

11. In subroutine C/)EFFI the imput cards concerning the profile 
losses of nozzle blades given in [6] are read, and the necessary 
interpolations are performed. 

12. In subroutine c/)EFF2 a similar process is followed for the 
profile losses of impulse blades. 

13 . In subroutine CALYA)6 the profile losses of nozzle blades are 
calculated with the coefficients computed in C/)EFFI . 

14. In subroutine CALYA2 the profile losses of impulse blades 
are calculated with the coefficients computed in c/)EFF2 . 

15 . In subroutine FIG6 the necessary interpolations for calcu- 
lating the secondary losses are performed. 

16 . In subroutine AINl/)S the overall loss coefficients with 
either method C 6J or 1 ^ 7 ] are calculated from blading geometry and 
flow angles. Subroutine AINL/)S must be called after subroutines 
C/)EFFI and Cj6eFF2, and calls subroutines CALYAj25, CALYA2 and FIG6. 
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The eight subroutines given below are concerned with method (ll) 
of the loss correlation. 

17. In subroutine TRAUPl the in^jut cards concerning the 
profile losses given by Traupel's method are read, and the necessary 
interpolations are performed. 

18. In subroutine XI^ the basic profile loss coefficients are 
computed with the coefficients calculated in TRAUPl. 

19. In subroutine CSIM Mach number effects on the loss 
coefficients are computed. 

20. In subroutine CID the effects due to the trailing edge 
thichness and nontwisted blades on the loss coefficients are con5>uted. 

21. In subroutine CSIW the end wall loss coefficients are 
computed . 

22. In subroutine CSIR the secondary loss coefficients are 
computed . 

23. In subroutine ALEAK the tip leakage efficiency losses are 
computed. 

2 k. In subroutine TRAUP2 the overall loss coefficients are 
confuted with Traupel's method CloJ from both blade geometry and flow 
characteristics. Subroutine TRAUP2 has to be called after TRAUPl, 
and calls subroutines XI^, CSIM, CID, CSIW, CSIR and ALEAK. 

25. In subroutine ZA1IX3S the overall loss coefficients are 
computed with the method of Lenherr and Carter L9] from the flow 
angles and the relative velocities. 

26. In subroutine BALJE the overall loss coefficients are 
computed with Balje's method [8j from blade geometry and flow 
conditions . 

e) Numerical Interpolation Subroutines 

27. In subroutine PARAB the coefficients of the parabolic 

interpolations used to approximate curves are calculated. 

28. In subroutine CHBFT the coefficients of the polinomial of 
arbitrary degree, which is the best fit of a given array of points 
in the Chebyschev sense, are calculated. 
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29. In subroutine YC the value of the polynomial determined 
in CHBFT, corresponding to a given value of the independent 
variable, is calculated. 

Description of Imput Data 

The imput data of the present con 5 )uter program consist of 63 + 

NSETS cards, NSETS being the number of conditions for which 

performance values are to be obtained. They can be divided in three 
sections. In the first, the constant values used in the methods 
of Ainley and Mathieson and of Traupel are introduced (these 36 
cards are read respectively in subroutines Cl^EFFI, Cj6EFF2 and TRAUPl). 

In the second, the specifications about the detail of output desired, 
the tolerances of the numerical solutions, the systems of correlaticn 
of the flow angles and the loss coefficients, the coefficients used in 
the equation of continuity, the gas properties and the turbine geometry 
are stated. In the third, the conditions for which performance values 
are to be found are introduced. 

The 36 input cards for the first section are given in Table I. 

They will remain unchanged for all turbines. Detailed description 
of the 27 cards composing the second section of the imput are given 
in Table II, where also a description of the third iir 5 >ut section can be 
found in cards #28 and # 29 . An example of the imput cards described 
in Table II is given in Table III. 

Description of the normal output 

In Table IV the output corresponding to the imput of Tables I 
and III is given. It can be seen that it consists of 1 + 3 NSETS 
sheets. In the first sheet the turbine geometrical dimensions and 
the choice among the different correlation methods are described. 

Then three sheets for every set of operating conditions are printed. 

The first of these sheets describes the flow conditions at station 1 
after the stator row; the second one the flow conditions after the 
rotor row; and the third one gives the main overall turbine characteristics 
and the mass averaged performance values. 
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Fig . I Coordinates, stations ond streamline nomenclature. 




Fig. 2 — Turbine thermodynamic process. 
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I CALCULATE INLET ROTOR CONDITIONS | 
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(AL0S2) 

( SEINE ) 



Fig. 4 Flow diagram 






Fig. 5 Velocity diagrams 



s 



AL = CAMBER LINE LENGTH 
C = CHORD 

0 = THROAT DIAMETER 

E = CURVATURE RADIUS 

S = SPACING 

T = MAXIMUM THICKNESS 

TE = TRAILING EDGE THICKNESS PROJECTED IN 
PERIPHERAL DIRECTION 

TN = TRAILING EDGE THICKNESS, NORMALTO FLOW 
DIRECTION 



Fig. 6 Blade geomety nomenclature 
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TABLE IV 



CORRELATION 

INDEX 

- 


POSSIBLE 

VALUES 




SPECIFICATIONS 




1 


A&M [6J or D&C [ 7 ] loss correlation 
method is adopted, depending on lAI 
value . 


IC0R 


2 


Balje [8] loss correlation method is 

adopted . j 

) 

! 




3 


j 

j 

L&C [ 9 ] loss correlation method is | 

adopted . 




4 


Traupel Cioj loss correlation method 


1 




is adopted. 


i 


0 

1 


If IC0R=1, AM4 [6] loss correlation 


lAI 


! 

1 


1 method is adopted. 

1 




1 

1 


t 

If IC0R=1, D&C [ 7 J loss correlation 
method is adopted. 




1 1 


i A&M L6] flow angles correlation is 

adopted . 
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Traupel LioJ flow angles correlation 


i 




is adopted. 
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Rotor is not shrouded. 
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Rotor is shrouded. 1 
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CORRELATION 

INDEX 


POSSIBLE 

VALUES 


SPECIFICATIONS 




0 


Loss coefficients are function of 
flow angles 


IINC 
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Loss coefficients are function of 
blade inlet angles, for negative 


‘ 




incidence 




1 


Loss coefficient C = Cq (see eq. (4?)) 


ICpL 


6 


Loss coefficient C = C (Y, pressure 
ratio) (see eq. ^l)) 




8 


Loss coefficient £ = g (see eq. (43)) 




1 


Loss coefficient in continuity equation 
C* = *5 CtoT 


ic<fm 


2 


Loss coefficient in continuity equation 
= ^PROF 




3 


Loss coefficient in continuity equation 

C* = 
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